Stochastic Pulse Switching in a Degenerate Resonant Optical Medium 
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Using the ideahzed integrable MaxweU-Bloch model, we describe random optical-pulse polariza- 
tion switching along an active optical medium in the A-configuration with disordered occupation 
numbers of its lower energy sub-level pair. The description combines complete integrability and 
stochastic dynamics. For the single-soliton pulse, we derive the statistics of the electric-field polar- 
ization ellipse at a given point along the medium in closed form. If the average initial population 
difference of the two lower sub-levels vanishes, we show that the pulse polarization will switch in- 
termittently between the two circular polarizations as it travels along the medium. If this difference 
does not vanish, the pulse will eventually forever remain in the circular polarization determined by 
which sub-level is more occupied on average. We also derive the exact expressions for the statistics 
of the polarization-switching dynamics, such as the probability distribution of the distance between 
two consecutive switches and the percentage of the distance along the medium the pulse spends in 
the elliptical polarization of a given orientation in the case of vanishing average initial population 
difference. We find that the latter distribution is given in terms of the well-known arcsine law. 
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I. INTRODUCTION 

Resonant interaction of light with active optical me- 
dia has given rise to one of the most fruitful areas of 
applied physics in having provided the foundation of 
numerous important physical effects over the past six 
decades p]-KL5|. and served as one of the basic mecha- 
nisms used in laser operation and optical amplifiers |16l - 
l2q |. While a fully quantum description of this interac- 
tion has also been developed [13, [21| , probably its most 
revealing description has been furnished by the semi- 
classical model provided by the MaxweU-Bloch system 
of equations [13, I2a - l25j . This model has helped to un- 
cover the fundamentals of the resonant interaction be- 
tween an electromagnetic field and a system of active 
atoms in the regime in which the field can be described 
classically and the medium by quantum mechanics, and 
in which a great number of the relevant experiments have 
been carried out [ll-[l5l|. In fact, many physical effects 
observed in these experiments, from the photon echo [2| 
and self-induced transparency [3, |j| to the chaotic laser 
dynamics |11| . have been explained using the Maxwcll- 
Bloch equations in the idealized two-level approximation, 
in which the light is assumed to be monochromatic and 
to interact resonantly with only two active atomic levels 
in the optical medium [13, l23l - |28j . 

In the simplest case of the two-level approximation, 
when the pulses interacting with the medium are much 
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shorter than the relaxation times of the medium, the 
MaxweU-Bloch system describing this interaction is com- 
pletely integrable [2^ • This feature was used to theoreti- 
cally explain three important phenomena: self-induced 
transparency, superfluorescence, and quantum amplifi- 
cation. The McCall-Hahn phenomenon of self-induced 
transparency [3[45[ — a medium whose atoms are initially 
in the ground state becoming transparent to optical 
pulses with the resonant carrier frequeuCT— was first an- 
alyzed using complete integrability in [29|, after many 
exact solutions hinting at possible integrab ility had been 
found in [30,[311- For superfluorescence J8l-[l0| — the gen- 
eration of optical pulses from the random fluctuations of 
the initial medium polarization in an excited medium — 
the linear stage was addressed in [84] , where the statistics 
of the delay time between the pumping of the medium 
and the pulse maximum were derived in terms of the 
statistics of the polarization fluctuations, and shown to 
be Gaussian. The fully nonlinear problem was subse- 
quently addressed using its integrable structure in ]33t - 
|35| , whose main result was the shape of the superfluores- 
cence pulse and its relation to the delay time. The fully 
nonlinear description of a quantum amplifier — incident- 
pulse amplification in an excited medium — was addressed 

in [H,!!!. 

An approximate description of the medium as having 
more than two levels, or degenerate levels, interacting 
with the light pulses propagating through it, is more 
physically realistic than the idealized two-level approx- 
imation. Such a description is desirable, for example, 
because effects such as self-induced transparency have 
also been measured for transitions between degenerate 



levels [38J. An important special active medium with a 
doubly-degenerate ground level and an excited level as 
its working levels is referred to as the A-configuration 
medium, so named because of the shape of the corre- 
sponding quantum transition diagram. The two types 
of atomic transitions in such a medium are stimulated 
by and emit light of two opposite circular polariza- 
tions [33, EO] ■ The complete intcgrability of the Maxwcll- 
Bloch systems describing light pulses interacting with 
this type of a medium was discovered in |4l| - t44| , and self- 
induced transparency was described. Superfluorescence 
and amplification of incident pulses via the transfer of 
energy from the initially excited medium to the pulse 
in A-configuration media were studied in [431 ^^^ f4a |- 
respectively. 

One feature distinguishing the A-configuration descrip- 
tion from the simpler two-level description is its ability to 
capture the polarization of the propagating pulses, and 
thus polarization switching which depends on the initial 
population of the two degenerate lower levels [43, |48| . 
Another distinguishing feature of the corresponding one- 
soliton solution is that it is a soliton only in the sense of 
being a potential in the direct scattering problem of the 
corresponding Lax operator that gives rise to a single- 
eigenvalue spectrum, but is not a solitary traveling wave. 
In fact, even in the simplest case of constant initial lower- 
level populations, its shape is only asymptotically sta- 
tionary. It has complex internal dynamics through which 
both its shape and velocity can change in space and time, 
and thus can refiect the light polarization switching. 

In this paper, using the corresponding integrable 
MaxwcU-Bloch system, we describe random polarization 
switching of pulses passing through a A-configuration 
medium induced by a disordered initial population of the 
degenerate lower sub-levels. The dependence of the prop- 
erties of an emerging light pulse on both intcgrability 
as well as randomness in the initial conditions appears 
already for the simpler ideal two-level optical medium 
through the phenomenon of superfluorescence [33i435| . 
Randomness in a two-level medium however appears to 
play a negligible role in self-induced transparency. Richer 
interactions between the integrable dynamics and ran- 
dom initial data emerge in a A configuration medium, 
as the ficxibility of populating the degenerate lower sub- 
levels makes it possible for self-induced transparency to 
take place in the presence of structural disorder arising 
from spatial fiuctuations in these populations. Such dis- 
order often results during the initial preparation of the 
resonant medium, and subsequently induces random po- 
larization switching of light pulses propagating through 
this medium. We will compute and analyze several sta- 
tistical properties of this nonlinear random polarization 
switching using exact results obtained with the inverse- 
scattering-transform technique for the A-configuration 
MaxwcU-Bloch equations. 

While, in reality, a pulse propagating through a reso- 
nant active medium encounters several sources of losses 
that make it decay on a number of relaxation scales. 



we have chosen its idealized lossless integrable Maxwell- 
Bloch description for two reasons. The first is that we 
aim to understand the fundamental features of the po- 
larization switching process in the course of this propa- 
gation, in particular how the pulse and medium param- 
eters affect its statistical properties. The second is that, 
at the current development level of the experimental in- 
strumentation, the situation in which the pulse width is 
much shorter than the relaxation times is achievable, so 
that our model should be realistic from the viewpoint 
of experiments. Therefore, we here consider the ideal- 
ized integrable model describing pulse interaction with a 
degenerate active optical medium in the A configuration 
with structural disorder introduced by an inhomogeneous 
distribution of the degenerate lower sub-level population 
in the medium. 

We present the polarization statistics for the one- 
soliton solution, both because we can obtain them ex- 
plicitly and because they yield a particularly transpar- 
ent description of the switching phenomenon. In our 
treatment, we use the classical polarization ellipse rep- 
resentation [43, [501, which has the advantage of being 
independent of time for the one-soliton solution; in other 
words, for the single-soliton pulse, the two angles deter- 
mining the shape of the polarization ellipse only depend 
on the location along the medium sample. We address 
the statistics of the pulse travel-time to a given loca- 
tion along the medium sample, the shape statistics of the 
polarization-ellipse at that location, and also the statis- 
tics of the switching between the left- and right-circular 
polarizations that a soliton pulse experiences while trav- 
eling along the sample. 

We explore the qualitatively different statistical dy- 
namical regimes that emerge for the cases when the initial 
degenerate-lower-level populations have (approximately) 
equal or unequal mean. In particular, we find that, when 
the lower levels are equally populated on average, the 
polarization lingers close to one of the two circular po- 
larizations for long distances but can forever switch in- 
termittently between the two with probability one. On 
the other hand, when the initial degenerate-lower-level 
populations along the optical medium have distinct aver- 
ages, the polarization after a few possible initial switches 
asymptotically approaches the circular polarization cor- 
responding to the transition between the on-average ini- 
tially less populated lower sub-level and the excited level, 
and no further switching occurs with probability one. 

The remainder of the paper is organized as follows. In 
Sec. ini we present the relevant background of the prob- 
lem. In particular, in Sec. Ill Al wc review the polarization 
ellipse representation of polarized light, in Sec. IIIBl we 
review the Maxell-Bloch equations that describe resonant 
interaction of pulses with a A-configuration degenerate 
active optical medium, and in Sec. Ill CI we we review the 
inverse-scattering-transform method and soliton solution 
used in the description of light-polarization dynamics. In 
Sec, mil we discuss the soliton statistics when the initial 
medium population is disordered in space, with the ap- 



proximate white-noise description presented in Sec. lIII Al 
and the statistics of the sohton travel time to a given 
point and the two angles determining the shape of the 
polarization ellipse at that point discussed in Sec. IIIIBI 
The statistical description of the polarization switching 
dynamics is given in Sec. IIV( and concluding remarks 
arc presented in Sec. [V] Appendix |X] further elucidates 
the appearance and role of the correlation length of the 
initial lower-level population difference along the optical 
medium. Appendix |B] contains a detailed calculation of 
two polarization-switching length statistics. 



II. BACKGROUND AND PROBLEM 
FORMULATION 

In this section, we describe the problem at hand 
by briefly reviewing the polarization ellipse description 
of polarized light, the three-level Maxwell-Bloch equa- 
tions that describe the propagation of monochromatic, 
elliptically-polarized light through a A-configuration ac- 
tive optical medium, and the soliton solution whose ran- 
dom polarization switching dynamics we will study in the 
rest of the paper. 



A. Optical Pulse Polarization 

The electric field polarization is among the light char- 
acteristics most sensitive to changes in the properties of 
the optical medium. This makes it a good potential tar- 
get for experimental investigation of the stochastic be- 
havior of the light pulses predicted in this paper. There- 
fore, in this section, we present a brief discussion of its 
main properties of importance to our subsequent discus- 
sion. 

Two well-established descriptions of light polarization 
are given in terms of the polarization ellipse and Poincarc 
sphere [4^,[53]- Here, we review the basic concepts of the 
polarization ellipse description, which we will use in the 
rest of the paper. (See also [43 ■) To this end, we con- 
sider a plane electromagnetic wave with frequency w and 
wave number fc, propagating in the positive x direction 
in the (a;,y, 2;,i)-laboratory coordinate frame, which has 
the form 



£{x,t) = Re|i?e*('' 

= Re { [SyGy 
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Here, Re denotes the real part of a complex number, Ey 
and Ez are the two components of the complex wave am- 
plitude, E, and e^ and e^ are the unit vectors in the 
(j/,z)-plane, perpendicular to the propagation direction 
of the wave. Defining the circular-polarization basis vec- 
tors e+ and e_ as 
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we rewrite the complex wave amplitude E in the circular- 
polarization components as 

E ~ EyBy + EzBz ~ E+ e+ + _E_ e_ 

= {\E+\e-'''' e+ + \E^\e"''e^) e'('^+'^), (3) 

where and 0-1-2-0 are the phases of the e+ and e_ 
electric field components, respectively. This yields the 
expression for the electric field 

- I F I 4- I F I 

£{x,t) = r= (cosipBy + simpez) 
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which clearly shows that this field traces out an ellipse, 
whose major and minor semi-axes have lengths (|i?+| -I- 
\E_\)/\/2 and (|i?+| — \E^\)/y/2, respectively, and whose 
major semi-axis subtends the angle ip with the y axis. 
The angle ip is called the polarization azimuth^ and takes 
on values < V' < tt. The ratio between the semi-axes 
of the ellipse is related to the ellipticity angle rj through 
the formula tanr/ = {\E+\ - |£;_|) / (|£:+| + \E^\). This 
angle takes on values — 7r/4 < rj < n/A; see Figure[T] The 
sign of rj represents the direction in which the electric field 
£{x,t) rotates along the perimeter of the ellipse. 

Special cases of Eq. (|4]) include linear polarization 
when T] = i\E+\ = |£^-|) and circular polarization 
when |?7| = 7r/4. In particular, the field £{x,t) is left- 
circularly polarized if 77 = — 7r/4, that is \E+\ = 0, and 
right-circularly polarized if 77 = 7r/4, that is |£'_| = 0. 

To compute the angles rj and tp, we proceed as follows. 
First, from Figure [1] we see that 
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Moreover, from ([3]), we find E'^E^ 
which implies 
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2iip 



tan 2-0 



E+E* 



E_E 



+ 



E+E*_ +E^E* 



(5b) 



Formulae ([5]) provide a complete characterization of the 
polarization ellipse. 

The concepts explained above hold equally well when 
the constant complex amplitude E is replaced by a com- 
plex amplitude E{x, t) that varies slowly compared to the 




FIG. 1. Polarization ellipse in the plane perpendicular to the 
direction of light propagation. 



plane carrier wave g^'^^""*-' . This is typically the case for 
the interaction of monochromatic light with a A configu- 
ration active optical medium [39J, |5l| , during which light 
pulses can be represented in the form 

f (.T, t) = Re { [£;+(.T, t) e+ + £;_ (x, i) e_] e'C'^^"*) } , 

(6) 
where e± are the circular polarization basis vectors ^ 
and E±(x,t) the complex envelopes of the two circu- 
lar polarization components of the light pulse inside the 
medium that vary slowly compared to the wavelength 
and time-period of the light. As explained in the next 
section, the two electric field polarization components 
interact with the two active atomic transitions in the 
A-configuration degenerate two-level medium. 



B. A-Configuration Maxwell-Bloch Equations 

Propagation of ultra-short, monochromatic, elliptically 
polarized light pulses interacting resonantly with a A- 
configuration, two-level, active optical medium, is, in 
the slowly-varying envelope approximation, described 
by the quasi-classical system of Maxwell-Bloch equa- 
tions [n, m, [m 

dH{t,x) dH(t,x) 



dt 



dp{t,x,X) 
dt 



dx 
[J,p{t,x,v)]g{v)dv, (7a) 



-XJ + H{t,x),p(t,x,X) . (7b) 



Here [•, •] denotes the matrix commutator. 

Equation (j7ap arises from the classical unidirectional 
Maxwell's equations for the electric-field envelopes with 
the displacement currents on the right-hand side. Equa- 
tion (|7bl) is the Liouville equation for the density matrix, 
describing in the present case a two-level quantum system 



with a doubly-degenerate ground level. The density ma- 
trix p, the reduced Hamiltonian H (without the diagonal 
part) describing the dipole interaction of the degenerate 
two-level system with the electric field, and the matrix J 
in Eqs. ^ are defined as 
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respectively. Here, E±{x,t) are the complex-valued en- 
velopes of the left- and right-circular polarization com- 
ponents of the light pulse as given in Eq. ©, p±{x,t, A) 
and /i(a:, t, A) the complex-valued medium-polarization 
envelopes, and n±{x,t,\) and Af{x,t,X) the real-valued 
population densities of the degenerate ground sub-levels 
and the excited level. The electric-field and medium- 
polarization envelopes, E±{x,t) and /5±(a;,i, A), arc as- 
sociated with the atomic transitions between each of the 
ground sub-levels and the excited level, while /i(x, t, A) 
is the contribution to the medium polarization by the 
two-photon transition between the two ground sub-levels. 
The parameter A describes the detuning of the atomic 
transition frequency from the exact resonance with the 
electric field, and g{X) is a nonnegative function with 
/_ dW dA = 1 which describes the shape of the spec- 
tral line due to the inhomogeneous broadening of the 
atomic transitions. The speed of light in Eq. (|7ap is 
non-dimensionalized to c = 1. In components, Eqs. ([7]) 
read 



dE± dE± 

H 

dt dx 

^-2^\p+ = \[E+{^f-n 
^-2iXp^ = \[E^{N-n 

^^l[E+*P-+E^p+*], 

dN 1 r 

— ==--[i?+p+ +E+ p+ 
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- E+p] , (9c) 
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(9e) 
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One condition for equations ^ (or ([9])) to be valid 
is that the pulse-width be much shorter than the time 
scale of the relaxation processes in the atomic system; 
as discussed below, in gases, the ratio between these 
two time scales typically ranges from 10~^ to 10"'^ jl7| . 
Also, as we already mentioned in the previous paragraph, 
Eq. (jTa)) describes unidirectional propagation. Poten- 
tial violation of unidirectionality is an important concern 
even for a non-degenerate two-level system, in which a 
spatially non-uniform density of active atoms can cause 



backscattcring of light. However, the most important fea- 
tures of resonant interaction between hght and two-level 
atomic systems are well-described within the unidirec- 
tional approximation [13, [H, H, [H-il [M [Slj. Since 
linear waves can be treated independently, bidirectional- 
ity must be taken into account only when the counter- 
propagating waves interact nonlinearly. Nonlinear inter- 
action, in turn, only becomes prominent when the wave 
amplitudes are sufficiently large and the characteristic 
time of the counter-propagating waves' overlap is longer 
than the characteristic onset time of the nonlinear inter- 
action. 

The amplitudes of the back-scattered waves are usu- 
ally small for two reasons: the low density of the active 
atoms, and the disorder in the lower sub- level populations 
of the A-configuration system leading to randomness in 
the phase of the back-scattered light. In particular, for 
the typical expected value of the electric dipolc, corre- 
sponding to the resonant atomic transition, of ~^ 1 De- 
bye, and the resonant transition frequency ~ 10^^ sec~^, 
the density of the active atoms < 10^^ cm~^ induces less 
than 2% of back-scattering according to the linear es- 
timates carried out in [S^l- (See also [5l[.) In addition, 
destructive summation of the back-scattered plane waves 
with random phases, which may result from the disorder 
in the lower sub-level populations of the A-configuration 
system, can also lead to an overall small amplitude of the 
back-scattered light. Finally, we should note that in prac- 
tical situations, the overlap time between two counter- 
propagating pulses is very short, due to the large value 
of the speed of light, and therefore so is the time of the 
nonlinear interaction between them. 

The density matrix, p, is Hcrmitian and its time- 
evolution can be represented by the formula p = UpoW , 
where U belongs to the group SU{3) and po is time- 
independent; this representation follows from Eq. (j7b[) . 
Thus, the three eigenvalues of the matrix p are con- 
served in time. Alternatively, we can find three inde- 
pendent conserved quantities for Ec^. (j7bp by computing 
the traces of the matrices p, fP' , and f(^ , three indepen- 
dent functions of the eigenvalues of p. Explicitly, these 
conserved quantities are: 



/i(a;,A) =7V + n+ +n_ = 1, 
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Note that unit normalization in Eq. (|10aP is chosen. 

The A-configuration Maxwell-Bloch system ([7]) con- 
tains two invariant two- level sub-systems [17J |. obtained 
by setting either _E-|_ = p^ — n+ = /z = or _E_ = 



p_ — Ti- = jjL = Q, which describe pure two-level transi- 
tions between the excited level and the — or -I- sub-levels, 
respectively. The light involved in either of these transi- 
tions forever remains circularly polarized. 



Polarization Dynamics in a A-Configuration 
Medium 



The solutions of the Maxwell-Bloch equations ([7]) can 
be obtained and analyzed via the inverse scattering trans- 
form starting with the zero-curvature representation [43| , 



dx 



V^ 
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m-'"' 



-iXJ + H 



(iXJ -H)^, 



ly-X 



(11a) 



g{v)dv ) $ 



(lib) 
where the symbol V stands for the Cauchy principal value 
of the integral, and the matrices H . p, and J are defined 
in formula ([S]). The 3x3 matrix $ is a simultaneous 
solution of both Eqs. piap and (jllbp . The compatibility 
condition of this system, Ux~Vt + [U, V] = 0, is equivalent 
to Eqs. d?]). 

The inverse-scattering transform is well suited to ad- 
dress the Cauchy problem for Eqs. ([7]) formulated along 
the entire real axis. Following [33, we thus introduce the 
asymptotically-mixed problem in which the Cauchy data 
represent a pulse incident at the point x = and defined 
along the entire <-axis. 



/OO 
\El{t)\dt <oo. 
-OO 



-OO < t < OO, 



(12) 



The asymptotic initial state of the optical medium is 
given at t = — OO by 



lim p±{x, t, A) ~ 0, 

i— f — OO 

lim N{x,t,\) =0, 



lim /i(a;, t, A) = 0, 

i— f — OO 



and 



lim n±{x,t,X)~— 1 ± a(a;, A) 

i— f — OO 2 



>0, 



(13a) 



(13b) 



with < X < L, where L is the non-dimensionalized 
length of the sample. Here, we have assumed that only 
the two degenerate lower levels are populated initially. 
The form of the asymptotic condition (|13bp follows from 
the normalization (|10a|) . which also implies that — 1 < 
a{x,\) < 1. 

In gases, the lifetime of the optical pulse ranges from 
IQ-^ to IQ-^ seconds, while the typical pulse-width is 
IQ-* seconds or shorter [13]. Therefore, the above ideal- 
ization of the initial time being at — c» is well justified. 



The initial conditions ((T2|) define the scattering prob- 
lem at the point x = for Eq. (|llap . which falls in the 
class of Manakov's scattering problems [53 • The evolu- 
tion of the scattering data in x can be obtained via Eq. 
(|llb[) . and the electric field envelopes E±{x,t) can then 
be recovered using a set of two Marchenko-type equa- 
tions [51|. The evolution equations for the scattering 
data corresponding to the most general asymptotic ini- 
tial state are listed in [4^ ; their derivation proceeds along 
the lines of the treatment given in [3a] for the single- 
polarization, two-level Maxwell-Bloch equations. 

As mentioned in the introduction, in [4l], using the 
inverse-scattering transform, a polarization switching 
mechanism was identified in the interaction of monochro- 
matic light with a A-configuration optical medium ini- 
tially satisfying the conditions P^ . In particular, as the 
pulse passes through an x-interval along which a{x, A) > 
is bounded below by a positive constant, the amplitude 
i?_(x, t) will grow in modulus towards a saturation value 
and E+(x,t) will decay, and vice versa if a(a;,A) < 0. 
If the initial state of the medium is prepared using un- 
polarized, incoherent light, the initial population differ- 
ence a{x, A) between the two lower sub-levels can be con- 
sidered a random function of x. Since therefore a{x^ A) 
changes sign in a random fashion, an optical pulse propa- 
gating in such a medium will experience random switch- 
ing of light polarization. Studying the statistical proper- 
ties of this random switching is the focus of this paper. 

Prom the inverse-scattering transform theory [5J, ISSf , 
it is well known that the asymptotic behavior of the so- 
lution will be determined by the discrete spectrum of the 
operator in Eq. (|llap — the iV-soliton solutions |48l. Isif. 
In fact, in the integrable Maxwell-Bloch-type equations, 
if the spectral line is not infinitely narrow (i.e., g{\) ^ 
(5(A), the Dirac delta function) the continuous radiation 
not only disperses away, but also becomes absorbed in the 
medium via Landau damping |29[ . If the discrete spec- 
trum of the incident pulse contains a single eigenvalue in 
the upper-half A-plane, i.e., Ai = "f + i/3, with /3 > 0, this 
pulse asymptotically reshapes itself into a single soliton 
{N=l). 

We address the case when the spectral width of the 
pump pulse is much broader than the width of the spec- 
tral line due to the inhomogeneous broadening, g{X). In 
this case, the initial populations can be considered homo- 
geneous within the width of the spectral line, and there- 
fore we can take 

a{x,X)=a{x). (14) 

The single-soliton solution is then given by the expression 

^±(x,t)=4i/3G±(x)e*®±("^*) 
X sech 



where the functions 
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control the amplitudes of the soliton components and the 
functions 

e± (x, t) === 27(i - x) + a[x± A{x)] - arg d± (15c) 

describe their phases; the real-valued coefficients a and 
T are given by 



5(^) 



\i — ly 



dv 



(16) 



for the given complex number Ai 
imaginary part, and 



7 + i/3 with a positive 



A{x) 



aiOdi, 



(17) 



describes the cumulative initial population difference 
a{x) along the medium sample up to any given position 

X. 

Equations (|15a|) and ()15bp show that both the electric- 
field components, E±{x,t), of the one-soliton solution 
consist of the same sech-profilc wave, with two differ- 
ent, x-dependent amplitudes, proportional to the func- 
tions G±{x) in Eq. (|15b[) . respectively. The maximal 
amplitude of each component equals 4/3. The tempo- 
ral width of the soliton equals 1/2/3. The constants d+ 
and d- determine the phase and position of the soliton. 
If the cumulative initial population difference diverges 
with increasing distance into the medium, A{x) — ^ ±oo 
as a; — > oo, one electric- field amplitude saturates while 
the other decays, which is the one-soliton case of the po- 
larization switching [43, |4^ . 

From Eq. (|T5)) . since g{v) > and Ai = j + i/3 with 
/3 > 0, we find for the coefficient t the inequality 



3(^) 



Ai-i^ 



rdj/ < 0. 



For the Lorentzian shape of the spectral line, 

7r(e^ + i^^) 
the coefficients in Eq. PB)) become 

7 /3 + e 



8 [72 + (/3 + ey 



i [72 + (/3 + e)2 



(18) 



(19) 



(20) 



Note that a = in Eq. (|16p if 7 = 0, which is also easily 
shown to be true for any even spectral line shape g{i^). 

The soliton speed and the phases of its components de- 
pend on the position x along the optical medium. Here, 
we compute the soliton speed as follows: At each posi- 
tion X, both components of the soliton reach their peak 



intensity at the time for which the argument of the sech- 
profile vanishes. This condition and Eq. (J15ap give the 
travel time of the sohton from when it is injected into the 
medium at x = until it reaches the point x as 



T{x) 



1 



Tx -\ — In 



2 \d+d. 



\d\ 



-In cosh ( 2tA{x) 



In 



\d- 



(21) 



From Eq. (|2ip . it is clear that the speed of the soliton 
thus satisfies the equation 

1 dr, 



-{x) = 1 — 

^'soiiton(a;) dx 2j3 



X tanh ( 2tA{x) + In ■ 



1 + a{x) 



(22) 



Using Eq. p^ and the fact that |a(a;)| < 1, we readily 
conclude that Vsoiiton(a;) < 1, i-e, that the soliton speed 
never exceeds the speed of light. If the initial popula- 
tion difference a is x- independent, the soliton velocity 
asymptotically behaves as 



"^soliton v^ 



1 



M(l-|a|) 
2/3 



To compute the polarization azimuth ^ and the angle 
of ellipticity 7] for the onc-soliton solution, we insert the 
components of the solution (|15ap into Eqs. ([5]) to obtain 

i: = -aA{x) + i arg {d*_d+) , (23a) 



(23b) 



sin 2-q — tanh ( 2tA{x) + In t- — - 



Note that these two angles are independent of the time t: 
at any point x along the medium, the light polarization 
remains constant in time as the soliton passes by. 

Note also that the polarization azimuth ip remains con- 
stant if the parameter a vanishes. Recalling that the 
eigenvalue Ai is the complex number Ai = 7 -|- i/3 and 
the remark after Eq. (|20|) . we see that this happens when 
7 = and so Ai = i/3, i.e., pure imaginary, provided that 
the spectral line shape g{X) is an even function. From 
Eqs. ([T5|) . its is easy to see that this case contains soli- 
tons with real-valued electric field components, which are 
obtained with the appropriate choice of the constants d± . 
In other words, when the spectral line shape g{X) is an 
even function, the polarization azimuth ip of all solitons 
with real-valued electric field components remains con- 
stant, and so is independent of the distance x into the 
medium. 



III. SOLITON DYNAMICS IN THE PRESENCE 

OF SPATIAL DISORDER IN THE MEDIUM 

POPULATION 

We now describe light propagation in the presence of 
spatial disorder in the initial population densities, char- 



acterized by the function a{x) in Eqs. (J13bp and (|14p . 
The spatial distribution of the initial population is de- 
termined by the manner in which the atomic system is 
prepared. In general, for the A-configuration with two 
degenerate levels, it is difficult to control the relative pop- 
ulations of the sub-levels during the preparation process. 
For example, if the system is prepared using unpolarized 
or partially polarized pump light, the relative distribu- 
tion of the sub- level populations will be random. We will 
be concerned with how this randomness induces random 
polarization switching in the one-soliton solution (|15ap . 



A. White Noise Approximation to the Initial 
Population Density Difference 

We assume the initial population density difference 
a{x) in the medium to be random and spatially homoge- 
neous in the statistical sense, and treat it as homogenous 
white noise with amplitude a superposed upon a mean 
(bias) b: 

(a(x)) = b, (24a) 

{[a{x) ~ b][a{x') - b]) = a^-d{x - x'), (24b) 

where (•) denotes ensemble averaging over the statistical 
ensemble of all possible realizations of the initial popula- 
tion difference a{x), and i5(-) is the Dirac delta function. 
The white-noise characterization ([M|) of the initial 
population density difference a{x) is consistent with the 
Maxwell-Bloch model ([7]) when the correlation length Lc 
of a{x) (discussed in more mathematical detail in Ap- 
pendix]^ satisfies three conditions. The first is that 



Lc > Ao, 



(25) 



where Aq is the wavelength of the light interacting reso- 
nantly with the transitions between the ground and ex- 
cited levels in the A-configuration medium under investi- 
gation. The second is is that Lc should be much shorter 
than the typical spatial pulse-width. 



1 



The third condition is 



ic< 



Lc < X, 



(26) 



(27) 



where x is the position of the observation point along the 
medium. 

As explained below, conditions (1^51) and (P5)) are neces- 
sary to make the modeling of the initial population differ- 
ence a(x) by random noise compatible with the Maxwell- 
Bloch equations ^ (or ((9))). Condition ((27)) is what 
allows the approximation of the true initial population 
difference a{x) by the idealized white-noise model (|24p . 
In fact, condition (P7| follows from condition ([^5)) in any 
realistic experimental device, which would be long com- 
pared to the soliton width 1/(3. We now proceed to dis- 
cuss the need and consequence of these three conditions 
in more detail. 



Condition ^5^ must hold because Eqs. © (or (O) de- 
scribe slowly-varying envelopes of the electric field and 
medium polarization components, and their carrier-wave 
oscillations are averaged out in the process of deriving 
these equations. The correlation length Lc must there- 
fore be sufficiently large compared to the wavelength Aq 
of the light interacting with the medium not to be aver- 
aged out as well. In other words, in order for the envelope 
approximation leading to Eqs. ^ to be valid simultane- 
ously with the assumption (|24bp , we must assume condi- 
tion (ps)) . 

Condition ^ should hold because Eqs. ^ (or ©) 
employ the approximation of unidirectionality. If the cor- 
relation length Lc of the medium with the initial popu- 
lation difference a{x) was comparable to or larger than 
the spatial width of the light pulse traveling through this 
medium, this random population difference could induce 
considerable backscattering of the pulse. Consequently, 
the pulse could be destroyed and the unidirectionality 
approximation violated. (Cf. the detailed discussion of a 
similar problem in [56|, |57[.) 

To understand the meaning of the condition (P7)) . let 
us recall that Eqs. p^ describing the soliton only in- 
volve the population difference a{x) through its cumula- 
tive spatial effect expressed by its spatial integral, A{x) 
in Eq. p7)) . In particular, when (P7)) holds, the integral 
A[x) can be well-approximated as 



A{x) = aW{x) + bx, 



(28) 



where W{x) is the usual Wiener process [S^. Recall 
that the Wiener process W{x) is for each x a mean-zero 
Gaussian random variable with variance x and probabil- 
ity density function 



Pwis;x) = 



1 



/2^ ^""^ V 2x 



(29) 



where s parametrizes the range of the random variable 
W{x). The representation (f28|) . which is equivalent to 
dMl), also applies to Eqs. dH]) and (gSl) for the soliton 
travel time and ellipticity angle, respectively. 

The parameter a in the approximation ([24| (or ([28]) ) 
is related to the correlation length Lc,and variance a^ 
of any given physical initial population difference a{x) 
through the equation a = \/2L^,aa, as follows from the 
discussion in Appendix 1X1 

We should remark, however, that the white noise ap- 
proximation (|24p does not make literal sense, as it vio- 
lates the constraint |a(x)| < I implied by the normaliza- 
tion ()10ap . More generally, it is not a valid description 
of physical quantities that depend on the initial popula- 
tion difference a{x) itself, such as the local soliton speed 
in Eq. (|22p . The precise description of such quantities 
would require a more detailed model of a{x), which we 
do not pursue here. However, state variables such as 
polarization variables and travel time, which involve the 
integrated effects of a(x), have their statistics well de- 
scribed by the white noise approximation (|24l) under the 
asymptotic conditions (P5|) through (|27p . 



In an experiment, the correlation length Lc of the pop- 
ulation difference a{x) would be approximately the same 
as the coherence length ic of the pump light used to pre- 
pare the optical medium. The related characteristic de- 
phasing time td = £c/c, where c is the speed of light (set 
to unity in our dimensionlcss coordinates), is determined 
by the width Ai'p of the spectral line of the light source 
as td ~ 1/Aup. Therefore, the coherence length ic is 
determined as 

ic ^ Ctd 



Ai 



Taking into account that the wavelength Ap is related to 
the frequency Vp as Ap = c/vp, we obtain 

dvp c 

dXp Xp 

and so, considering just the magnitudes of Avp and AAp, 
we find that 

Avp = -72 AAp. 
Ap 

From this formula, it finally follows that the coherence 
length ic is determined as 

X-n 

fl ^ . P 



AXr, 



Here, Ap is the average wavelength of the pump light and 
AAp is the characteristic width of the pump light-source 
spectral line (in terms of wavelength) . 

To demonstrate experimental feasibility, we recall that 
for Ti-sapphire lasers, for example, these parameters are 
AAp ^ 5nm and Ap ^ SOOnm, therefore ic ~ lO^nm = 
0.1mm [59|. For a typical A-configuration transition pair 
in the visible regime (e.g., in sodium vapor, the wave- 
length corresponding to the transition is Aq ~ 600nm), 
this argument shows that Xq <S^ ic ''^ Lc on the one hand, 
and that a several-centimeters long experimental device 
is clearly sufficiently long to capture the desired statisti- 
cal effects. Therefore, both conditions (PS|) and (P7)) can 
be satisfied simultaneously in this case. 

Noting that the soliton travel time T{x) and the ellip- 
ticity angle ri{x) depend on the initial population differ- 
ence only through the product tA(x), defined in Eqs. (|I6p 
and (jI7p . respectively, we here identify three fundamen- 
tal length scales associated with the dynamics of this 
quantity, and thus, through ri(x), also the polarization 
switching. First, as seen from Eq. 



1 



(30) 



is the length scale over which the deterministic bias b in 
Eq. (j24ap induces a significant change in tA{x). Second, 



L, 



(31) 



is the length scale over which the random component 
of the initial population difference fluctuations in the 
medium, given approximately as aW{x) in Eq. ([28| . cre- 
ates a signficant change in tA{x). Finally, 



^flu 



7^2 

La 



a 
6^' 



(32) 



is the length scale before which random fluctuations dom- 
inate the effects of the deterministic bias, and after which 
the opposite is true. Because of Eq. ([5^ . these three 
length scales must obey one of the two following order- 
ings: 

La < Lb < Lfluc or Lfiuc < Lb < La- (33) 

Note that 



b ^ oo. 


Lfluc -^ oo 


as 


6^0, 


(34a) 


Q -> OO, 


ifluc -> 


as 


a^O. 


(34b) 



Note also that the lengths La and Lb depend on the initial 
medium parameters a and b in Eqs. (|24p as well as the 
soliton parameter r in Eq. (|16p . while ifluc only depends 
on a and b. 

The polarization azimuth 'ip{x) depends on the disorder 
of the initial medium occupation numbers through the 
product aA{x) which, as we will see in Scc. lIIIB 2al does 
not require length scales analogous to La and Lb- 

The question that we need to answer is how the polar- 
ization of the pulse behaves at large distances x into the 
medium. If the initial difference a{x) between the pop- 
ulations of the degenerate ground states of the medium 
exhibits a bias, {o:{x)) = b, it is reasonable to expect 
that the soliton will eventually evolve into a single cir- 
cular polarization, which will depend on the sign of this 
bias. If no such bias exists, that is, {a{x)) — b — 0^ then 
it is reasonable to expect that the soliton will switch in- 
termittently between left and right circular polarizations 
over large distances. In the forthcoming sections, we will 
confirm this intuition explicitly. 
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FIG. 2. The two circular-polarization amplitudes, E+ and 
E-, of a real- valued, polarized- light soliton pulse propa- 
gating through a randomly-prepared A-configuration optical 
medium, as described by Eq. (|15al) . The model of the ini- 
tial population difference a{x) was sampled from the prop- 
erly rescaled beta distribution with mean b = {a{x)) — 0.001, 
variance a% — 0.71, and coherence length Lc — 1.8. The 
width of the Lorentzian spectral line in Eq. (|19[) is e = 0.1. 
The soliton parameters are /3 = 1/3, 'y = 0, d+ — d- — i. 



Soliton Travel Time 



B. Soliton Statistics at Fixed Observation Point 

In this section, we calculate the statistics of the soli- 
ton travel time T{x) in Eq. (PT|) and the two angles that 
determine the dynamics of the polarization ellipse, i.e., 
the polarization azimuth ip and angle of ellipticity rj in 
Eqs. (j23ap and (j23bl) . Throughout the section, we em- 
ploy the white noise approximation (j24p of the initial 
population density difference a{x) (or, equivalently, the 
Wiener process approximation (^5)) for its spatial integral 
A{x)). From the previous section, we recall that this re- 
quires the observation point x to be sufficiently far into 
the medium in comparison with the correlation length Lc 
of the function a{x) (cf. Eq. (gT])). 



As we recall from Section fllCi the soliton travel time 
T{x), given by Eq. (PT|) . is the time needed for the 
peak of the soliton to reach the observation position x. 
The time T{x) is a random function with the random- 
ness arising solely from the integral A{x) of the initial 
population density difference a{x), defined in Eqs. (jI7p . 
(J13b[) . and (fTi|) . respectively. As in Section fill Al we as- 
sume the Wiener process representation ([28)) for A{x), 
i.e., A{x) = aW{x) + b. In the calculations below, for 
each fixed x, we parameterize the range of the random 
variable W{x) by the variable s. 

a. Mean and Variance of the Soliton Travel Time 
The mean and variance of the soliton travel time T{x) 
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can be expressed as the integrals 

r(x, s) 



(r(x)) 



\/27r.T 
T{x,s 



exp 



ds, 



(35a) 



{T{x)) 



'2tix 



■ exp 



-2^ ' ^^- 



(35b) 



where T(a;, s) is defined as ^(a;) in Eq. (PT|) with A(a;) 
replaced by as + 6a;. 

In general, the integrals in Eqs. ((35|) can only be eval- 
uated numerically. For sufficiently large distances, how- 
ever, we can exploit the formula 



Incosh(u) 



ln2 + 0(e-^l"l) 



(36) 



valid for \u\ '3> 1, in Eq. ((2T|) . to approximate T{x,s) in 



Eqs. (|35j) as 



Tix,s) 



1 



Tx -i — In ■ 



\d+d^ 



\d. 



2r{as + bx) + In — — - 

Qj — 



+0(e- 



-\2T(as+bx)\\ 



(37) 



In this and the following asymptotic statements, we will 
treat ln(|(i+|/|d_|) as a fixed quantity of order unity 
(which, in connection with Eq. (j23b[) . means the pulse 
polarization at the entrance to the medium is not close 
to circular); otherwise, additional length scales involving 
ln(|d+|/|d_|) would appear in the error estimates. The 
mean soliton travel time deep into the medium can then 
be expressed as 



(Tix)) 



1 

2^ 



Tx{l-\b\) 



2^"M.P 



\d-\ 



\d-\ 



+0 (e-^/-^^ + (a;/L6),e-^/-^"-) , 



for X > ib,iflu 



/2" I ^^1, \d+d 



-O ((x/L„)-l/2 + (x/L„)l/2(a;/Lfl„,)) , 



2/3 



\d- 



(38a) 



for La^ X <€. ifiu 



and its variance as 



al{x) 



+ O (e"2-/i^ + (a;/Lb)2e-"/^«-) , for x » Lf,, Lfl„c, 

\d+\ 



4/32 

{tx - 2)t'^ a^ X 4-7r 

47r^2 167r/32 V'"" M- 

^+0 ((.T/La)-'/' + (.T/Lj(x/Lfl,e)^/') 



In- 



(38b) 



for La <C a; <C Lflu 



where the parameter t is defined in Eq. (|16p. 

Note that since r < due to Eq. ([T8| . the expecta- 
tion value of the soliton travel time T{x) in Eq. (|38ap 
increases linearly for large observation-point distance x. 
Recalling from Eqs. (|30)) . (|3T|) . and ([32)) that the regime 
a: ^ Lfluc is dominated by the random fluctuations in 
the initial population difference, whereas the drift due 
to the bias b in the difference dominates at length scales 
X 3> Lflucj we note that the first line in each display cor- 
responds to the drift-dominated case while the second 
line corresponds to the noise-dominated case. In partic- 
ular, the case x ^ Lh,ifluc only corresponds to nonzero 
average initial population density difference b = (a(a;)), 
while the case L^ <C a; <C Lfluc also contains the case 
6 = 0. Because of the length-scale relationships (|55)) . the 
regimes considered provide a comprehensive description 



I 

for the soliton travel time deep in the medium. 

b. Probability Distribution of the Soliton Travel Time 
The cumulative distribution function 



Fr{t;x) = Prob|r(a:) < t\. 



(39) 



for soliton travel time, T{x) in Eq. (HI]), from the en- 
trance of the medium to a given position x, can be 
computed in the Wiener process approximation using 



Eqs. (1211) and (p8|) . which yield the expression 

Fr{t;x) =Prob|l^(a;) < u^{x,t) or W{x) > u+{x,t)\ 
l,Prob|Ty(a;) < U-{x,t)\ 
+ Prob|vF(x) > u+{x,t)\ 

where 



u± (x, t) = - — ( — 26t.t — In ——7 q= cosh ^ ^(x, t) 
2ar V \d-\ 



and 



(40) 



k[x, t) = exp ( -4/3(t ~x)- 2tx - In ,_, + ' 



M+P + M- 



(41) 
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Since W{x) is a normahy distributed random variable 
with mean and variance x, we find 



1 ^ fu^(x,t) 
1 + - erf y=^ 

2 V V2^ 



Fr(t;x) 



ierff^^M 
2 V V2^ 



for < i < imax(a;), 

1, for t > tmax(a;), 

(42) 
where tmax(a;) is the upper bound on the soliton travel 
time to a position x, given by 

, , f T \ 1 , 2\d+d_\ 

W(a:)=(^l-^jx--ln — 



and the error function is defined as 



2 f^ 

V '"' Jo 



(43) 



The probability density function PTit] x) of the soliton 
travel time T{x) is given by the partial derivative of the 
cumulative distribution function F-Y{t; x) with respect to 
t. After differentiation of Eq. (|42|) and some algebra, we 
arrive at the expression 



pr{t;x) 



2 l3K{x,t) 

T^x a\T\.jK'^{x,t) - 1 



exp 



.0, 



u\{x,t) 
2x 



exp 



u'^_{x,t) 
2x 



for < f < tmax(a;), 



otherwise. 



(44) 



with u±(x, t) and K(a;, i) as in Eqs. ()40p and ()4ip . respec- 
tively. 

A sample plot of the probability density function 
PT{t'i x) of the soliton travel time T{x) is presented in 
Figure [3l Note the local maximum emerging from the 
boundary at i = imax(a;) at large values of x, i.e., loca- 
tions deep into the medium sample. 



I 

sian form 

Pip{s;x) = 

X exp ■ 



1 



/2ttx aa 



s — i a.Tg{d*_d+) + abx 



2xa^a^X) 



with mean 



(45) 



2. Polarization Variables 

To describe the statistics of the polarization azimuth 
ip and the angle of ellipticity rj, we again use the Wiener 
process approximation (|28p and replace the function A{x) 
in the expressions (j23p for the polarization variables with 
aW{x) + bx. The statistics can then be obtained as fol- 
lows. 

a. Polarization Azimuth Statistics As the polariza- 
tion azimuth ip{x) is a linear function of the Wiener pro- 
cess W{x), it itself behaves like a Brownian motion with 
drift —ah and diffusion coefficient ^a'^a'^. That is, its 
probability density at any position x is given by a Gaus- 



(■0(2;)) = -abx + - aTg{d*_d+) 



and variance 



(7l{x) = a^a^x. 



Note that when ct = 0, 



(46) 



(47) 



p^(s; a;) = (5 ( s - - arg(dld+) 



where (5(-) is the Dirac delta function, i.e., the dynam- 
ics of the polarization azimuth becomes constant, as was 
mentioned at the end of Section III Al In particular, the 
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FIG. 3. Probability density function, pr{t',x), with /3 = 1, 
7 = 1, e = 0.1, d+ = i, d- = i/3, a = ?>, h — 0.8. 



CiH 




100 



s^ 




100 



FIG. 4. Probability density function, p-^{s;x), with /3 = 1, 
7 = 1, e = 0.1, d+ = i, d^ — i/3, a = 1. Top: no bias, 6 = 0. 
Bottom: nonzero bias, b = 0.3. 



polarization azimuth is constant for all solitons whose 
electric-field envelopes are real- valued, so that for such 
solitons, the polarization ellipse does not rotate. 

b. Mean and Variance of the Ellipticity Angle For 
the angle of ellipticity t]{x), just as for the soliton travel 
time T{x), the initial medium population difference aix) 



again contributes to the evolution of rj{x) in the x direc- 
tion through the Wiener process tA{x). Therefore, the 
length scales given by Eqs. ([5n| . ([5T|) . and ([5^ are again 
relevant here. 

Taking H(a;, s) to be defined as the expression for sin 27] 
in Ec^. (|23bp with A{x) replaced by as + bx, i.e., 



(48) 



H(.T, s) = tanh 2T(as + bx) + In '—-^ 

\d\ 



the expectation and variance of sin 2?; can be expressed 



as 



(sin277) = (H(x,M^(x))) 



H(a:,s) /^ s'\ , 

— , exp — -— as 

(49a) 



'''sin 2j; 



(H(a;,.s)-(sin(2?7)))2 



■ exp I as. 

/2^ V 2x7 

(49b) 
Although the integrals in (j49| must, in general, be 
computed numerically, they can be evaluated asymp- 
totically at locations deep into the medium. When 
La <^ X -^ .^fiuc, rescaling the integration variable as 
s = y/xs' , we see that 



H(a;, y/xs) ~ tanh 2T{a\/xs + bx) + In -— — r 
~ sgns' 



(50) 



while the Gaussian integration factor becomes indepen- 
dent of X. We can therefore evaluate the La -^x -^ ifiuc 
asymptotics of both the mean and variance by the asymp- 
totic replacement ((50)) . and find that (sin 27; (x)) ~ and 
^lin2,M) ~ 1 for Lq < a; < Lfluc- 

On the other hand, if x ^ Lb, Lfluc, we have under the 
same rescaling 



H(x, y/xs') = tanh ( 2T(as/xs' + bx) + In [-i| 



sgn6. 



sgn6 and a^-^ 



and so we find that (sin 2?/) ' 
large x. 

To recapitulate, we have that for large x 



sin 2r] 



(sin 2?7) 



sin 27/ 




-sgn5, x>Lb,Lfiuc, 

0, La'^X <t: ifluc, 



X » Lb,Lfluc, 

La < X < Lfl„ 



(51) 
for 



(52a) 



(52b) 



Again, we must recall from Eqs. ([50)1 . ([5l]) . and ([5^ that 
the case x 3> Lf,, Lfluc corresponds to only nonzero av- 
erage initial population density difference b = (q;(x)), 
and the case L^ <C x ^ Lfluc also contains the case 
6 = 0. We therefore see that all solitons with nonzero bias 
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b = {a{x)) will, with probability one, collapse into a per- 
manent circular polarization. When b = 0, the large x re- 
sults for the mean and variance of sin 2r] in Eqs. (|52|) show 
that ([sin 277(0;)]^) = cr'^ia2rii^) + {sinrj{x))'^ —)■ 1, which 
along with (sin277) — > and the inequality | sin2?7(a;)p < 
1, implies that sin2r/(a;) must converge at large x to a 
random variable concentrated at ±1. with equal weight. 
In other words, for large x, the light polarization becomes 
one of the two circular polarizations with equal probabil- 
ity. 

c. Probability Distribution of the Angle of Ellipticity 
The large x asymptotic behavior of the ellipticity angle 
77, described at the end of the previous section, can also 
be seen from developing the exact formula for the cumu- 
lative distribution function F^(s;a;) of 77, which can be 
computed from the expression for 77 in Eq. (j23bp via a 
formula analogous to Eq. ([M)) . In this way, we find 

F^{s;x) = Probi W{x) > - — ( tanh^^(sin2s) 



-2Tbx — In ■ 



\d- 



(53) 



Since W{x) is normally distributed with mean and vari- 
ance a;, we compute 



Fvi^'X) ^2^^^^ 



X, tanh (sin2.s) 

' 2aT V 



-2r6x — In 



\d-\ 



(54) 
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FIG. 5. Probability density function, prj{s;x), with P = I, 
^ — 1, e — 0.1, d+ — i, d- — 3i/4, a — 1. Top: no bias, 
6 = 0. Bottom: nonzero bias, b — 0.3. 



where x(a;, u) = erf (u/v2a;), with the function erf(-) 
defined in Eq. (gS]). After differentiating Eq. dM} with 
respect to the parameter s, we derive the probability den- 
sity function 

1 

Pr,{s;x) 



12-Kx alrl cos 2s 



X exp < 



tanh ^(sin2s) -2T6a;-ln|rf+|/|d_ 



2 ^ 



SaV^x 



> . 



(55) 



Note that formulae ([54|) and ((55)) arc only valid for 
— 7r/4 < s < 7r/4, the range in which the ellipticity angle 
r\ is defined. 

Very far into the active medium, that is, for large val- 
ues of X, the distribution function p,,(s;x) in (|55p clearly 
attains very small values at all s away from s — ±t:/A, 
while at s = ±7r/4, it exhibits singularities. Indeed, as 
remarked in Section UlIB 2 bl the large x asymptotics of 
the mean and variance of the angle of ellipticity 77 (x) im- 
ply that its probability distribution must concentrate at 
one or both of the values s = ±7r/4 corresponding to the 



two circular polarizations: 



p,,(s;x) 



1 r 
2 



S s 



+ Sis + 



S 
6{s 



TT 

^+4 



, 6 = 0, 
6>0, (56) 
6<0, 



where 6{-) is the Dirac delta function. 

The discussion in the preceding paragraph shows that 
for very large distances x into the medium, the soliton 
will mostly be confined to one of the two circular polar- 
izations. For nonzero average initial population density 
difference, {a{x)) = 6 7^ 0, this polarization is fixed by 
the sign of (a(x)) = b and with probability one even- 
tually stops switching. For (a(x)) = 6 = 0, at large 
distances into the medium, the soliton stays in one of the 
two circular polarizations for most of the time, switch- 
ing intermittently between them. The dynamics of the 
switching will be discussed in Section HVl 



14 



IV. DYNAMICS OF POLARIZATION 
SWITCHING 

Having developed explicit forniulas for the soliton 
statistics as functions of depth into the optical medium 
when the initial population difference a{x) is random, 
we now provide a brief quantitative description of the 
dynamics of polarization switching. We begin in Sec- 
tion IIVAI by identifying some key length scales to de- 
scribe the essential features of the polarization dynamics. 
Then, in Section FlVB( we present some analytical results 
for the polarization switching dynamics in the Wiener 
process approximation. 



A. Length Scales of Polarization Switching 
Dynamics 

Because the asymptotic states of light pulses inter- 
acting with a A-configuration medium are given by the 
two circular polarizations, crossing the linear polarization 
represents a key reference point on the pulse trajectories. 
In particular, the key elementary stages of the polariza- 
tion switching can be cast in terms of this crossing: The 
pulse will generally evolve from a linear (or elliptical) 
polarization to a nearly circular polarization, and then 
possibly eventually return to a linear polarization, from 
which it could return to its previous or the opposite cir- 
cular polarization. The characteristic length scales corre- 
sponding to the transitions between a linear and circular 
polarization and the successive returns to a linear po- 
larization need not be the same, because the pulse may 
reside near a circular polarization for long distances be- 
fore returning to a linear polarization. Moreover, for po- 
larization switching in the presence of a nonzero bias, 
b = {a{x)) 7^ 0, a pulse will eventually remain in one 
circular polarization forever, and so it is important to in- 
troduce a characteristic distance after which there is no 
more switching. 

In view of the discussion in the previous paragraph, 
we can identify three distinct lengths of interest associ- 
ated with the light pulse polarization switching process: 
the distance between successive switches, the length scale 
over which the switching process manifests itself when it 
does occur, and the distance into the medium over which 
switching continues. One can distinguish these length 
scales more precisely by defining the following random 
distances: the switching transition distance Xt,.a over 
which the light pulse polarization evolves from a linear 
state (r/ = 0) to a nearly circular polarization state of ei- 
ther orientation (|?y| = ir/A — rjc), the interswitch distance 
Xint over which the light pulse polarization evolves from 
a linear state (77 = 0) to a nearly circularly polarized state 
(1 77 1 > 7r/4 — i]c for some fixed rjc > 0) of either orienta- 
tion and back to a linear state {rj = 0), and the switching 
region depth Xdop beyond which the light pulse polar- 
ization remains forever in one of the circularly polarized 
states I77I > 7r/4 — r]c for all greater distances. As noted 



above, the interswitch distance Xint is not necessarily the 
same order of magnitude of the switching transition dis- 
tance Xtra because X-mt also includes the distance over 
which the soliton remains in a circular polarization before 
returning to a linearly polarized state. 



B. 



Polarization Dynamics in Wiener Process 
Approximation 



From the equation (|23bp describing the spatial depen- 
dence of the ellipticity angle r] on the position x along the 
medium sample, one can see that the distances Xtm, -'^int, 
and Xdcp depend solely on the level-crossing properties of 
the Wiener process 2rA(a;)+ln(|(i+|/|(i_|). In particular, 
computing Xtia, ^intj and X^^p is equivalent to finding 
the positions along the medium sample for which this 
process first reaches the absolute value tanh~^ [cos(277c)] 
after having passed through the origin, first returns to 
the origin after such an excursion, and remains further 
from the origin than this absolute value for all subsequent 
X. We consider separately the case in which the initial 
population density difference a{x) in the medium has no 
bias (b = {ct{x)) ~ 0) and in which it does have bias 
{b={a{x))^0). 



1. Case of no Medium Bias 

When the initial population-density bias vanishes, 
Eq. (j52p implies that the polarization observed at any 
given position deep into the medium is likely to be circu- 
lar, with probability 1/2 for each orientation. From a dy- 
namical perspective, in fact the polarization switches in- 
finitely often, arbitrarily far into the medium, with prob- 
ability one, because the Wiener process is recurrent in 
one dimension. Consequently X^cp = 00 with probability 
one. The polarization does indeed reside over great dis- 
tances within one or the other circular polarization state, 
punctuated occasionally (but persistently) by switches 
(over relatively short distances) into the opposite polar- 
ization. 

Finding the statistics of the distances Xtra and Xjnt is 
equivalent to finding the corresponding distance statistics 
for the Wiener process 2TaW{x) without drift. There- 
fore, we can apply the well-known first-passage-time 
formulas (see, for example. Sec. 7.3 in [60| or Sec. 
2.8 in [6l|), including their convolution computed via 
Laplace transforms, as described in Appendix [B] to ob- 
tain the following expressions for the probability den- 
sity function pxtra {^) for the transition switching dis- 
tance Xtra, and the probability density function px^^^ix) 
for the interswitch distance Xjnt, using the definitions in 
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Sec.|Tv3 



PXt„ (x) 



2itr 



y^ (4n + 1) exp 



■-(4. + !)^^ 



, (57a) 



PXi„t (a;) 




-(4n + 2)^^' 



16n 



2 -^tra 



2a; 



where 



Li, 



2a|r 



taiih ^ (cos(277c)) 



(57b) 



(58) 



set the corresponding switching length scales (cf. the 
distance La in Eq. (pij) .') These two distributions are 
depicted in Fig. |6l 
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FIG. 6. The distributions pxtr!i.{^) (solid line) and pXi^^ioo) 
(dashed line) for a medium with no bias, 6 = {a{x)) — 0, 
computed using Eqs. (|57p with Ltra = 4.53. Note the faster 
decay of the distribution pxtis_{x)- Inset: The distributions 
pxtra (x) (solid line) and Px^^^ (x) (dashed line) for a medium 
with strong bias, computed using Eqs. (|63|l and (|66p . respec- 
tively, with a = 1, b — —0.5, and Ltra = 4.26. 



Note that the interswitch length and transition switch- 
ing length appearing in the probability density functions 
are the same, but the probability distributions are quite 
different. In particular, the transition switching distance 
has finite mean and variance: 







FIG. 7. Polarization state in a medium without bias under the 
Wiener Process Approximation (a = 1.00, fe = 0). The length 
scales for the transitions and distances between switches are 
here Ltra ~ Lint = 4.53, but note that the distance between 
switches has higher probability to take values large compared 
to this typical length scale. 



On the other hand, the probability density function for 
the interswitch distance Xi^t is so slowly decaying that 
Xjiit has an infinite mean, even though we have identi- 
fied a finite interswitch length scale ([55)1 . The meaning of 
this is that while many polarization switching events do 
take place with interswitch distances comparable to Lint, 
on occasion a much longer distance is observed between 
polarization switches, and these rare events still have a 
large enough probability to imply an infinite mean inter- 
switch distance. The transition switching distances are 
however much more likely to be on the order of Ltra, as 
can be seen from Eq. ((59|). That is, the polarization will 
fairly often tarry in a circular polarization state for a dis- 
tance much larger than Ljnt, but it will usually move out 
of the linear polarization state over the length scale Ltra, 
as shown in Fig. [T] This is naturally reflected in the po- 
larization statistics developed in Eq. (|52p , which indicate 
that the polarization will, deep in the medium, tend to 
be in one of the two circular polarization states. 

Another quantifiable statistic for the case of no bias 
in the initial population density is the fraction $ of the 
length over which the polarization ry takes a certain sign. 
One can show (see [63 or Sec. 1.4.4 in [63) that its prob- 
ability distribution is given by the arc sine law, which has 
the feature of having a rather large probability to take 
values near </> = or = 1 , meaning that even though the 
medium is unbiased, an individual realization is rather 
likely to spend most of its observed time in one or the 
other circular polarization state: 



Prob(<l' < 



(^tra 



Lt 



-L 






(59) 



p^{(t>) 



3in ^ V^ : 



7rV'/'(l -</-)' 



P*( 
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2. Case of Nonzero Medium Bias 

If the initial population of the atomic ground levels in 
the medium does have a bias (6 = (a(x)) ^ 0), then po- 
larization switches due to random fluctuations can occur 
for a while, but the bias in the population density of the 
optical medium will with probability one eventually col- 
lapse forever into the preferred circular polarization state 
associated with the sign of the initial population density 
bias h = (cx{x)) [48|. More precisely, if the bias b has 
opposite sign to that of the eUipticity angle. 



1 
?7o = - sm 



tanh In 



Mi 

\d- 



(60) 



of the pulse upon entering the medium, the polarization 
of the soliton will proceed through the following stages: 
achieving a linear polarization for the first time after a 
distance Xnn^Oj rnoving into a favored circular polariza- 
tion state over a subsequent distance Xtia, and never 
leaving this ultimate circular polarization state after a 
subsequent distance Xfluc- 

The polarization switching depth Xjop would then be 
the sum of the lengths corresponding to these stages 
^dop = ^iin,o + ^tra + ^fluc- The iutcrswitch distance 
Xint is not well defined for the case of a medium with 
bias because eventually the soliton will stop switching. 
If the bias were of the same sign as the initial polariza- 
tion state, then Xn^fi would be meaningless and should 
just be treated as in what follows, and the probability 
distribution for Xtra would be changed only by replacing 
the length scale for Ltra in Eq. (|M1) below by 



J-itra. — 



1 



2|t6| 



tanh (cos(2?7c)) 



In 






First we develop formulas for the statistics of these 
lengths, then discuss the qualitative differences between 
soliton polarization evolution in a medium with weak bias 
and with strong bias. The probability density function 
PXun i^) fo^ the distance Xn^fi until a linear polarization 
?7 = is first reached can be expressed through another 
first passage time formula for Brownian motion with drift 
(see Sec. 7.5 in je^l), 



Px„„,o(a;) 



\b\L 



lin,0 



where 



and 



i\/2t 



^lin.O 



; exp 



fe^(iiin,o-a;)^ 



2a^x 



1 ,Jd,\ 



2\Tb\ \d^\ 



(-'^Iin,o) — L\in,0, 
2 O, Llin,0 



'Xu 



62 



(61) 

(62a) 
(62b) 



After reaching the linear polarization state, the polar- 
ization will tend to move toward its favored circular po- 
larization state, reaching it after a further distance Xtm 
which has probability density function (see Sec. 7.5.5 
in [60|) 



/ N Ml^tra 

Pjftra [X) = 7^^ exp 



a\/2Trx^ 



b\Lt,^-xf 



2a^x 



where 



and 



-^tra = -T—r: tanh ^[cOs(2?7c)] 

2\Tb\ 



{Xt, 



L 



-k. 



tia; 
J^tra 



(63) 

(64) 

(65a) 
(65b) 



Once the polarization achieves a value \ri\ > tt/A — 7]^ 
corresponding to a circular polarization with the same 
sign as b, it will continue to revisit the linear polarization 
state r/ = over a further distance Xf\ „c., fo r which a last 
passage time formula (see Sec. IV. 5 in [6J|) shows that it 
is distributed as LfiucG^, where G is a standard Gaussian 
random variable with mean zero and unit variance, and 
Lfiuc = o? Ib"^, as defined in Eq. ([32]) . sets the length scale 
over which switching behavior continues. In other words, 
^fluc is governed by a x^ distribution with one degree of 
freedom, and has probability density |62| 



PXnuAx) 



2'KX 



exp 



Lfiv 



(66) 



and mean (Xa-nc) = -^fluc- The distributions PXtr^i^) in 
Eq. ([63)1 and pxnua i^) ^^^ depicted in the inset of Fig. |6l 
The qualitative character of the soliton trajectory will 
depend on the ratio 



(-'^flu 



(^tra 



(67) 



where, from (HS I, {X ruc) = ifiuc = a'^/b'^, and (Xtra) = 
Ltra is given in (|65ap . The numerator describes the dis- 
tance over which the polarization continues to fluctuate 
between the two circular polarizations, whereas the de- 
nominator characterizes a single transition from linear to 
circular polarization. From the formulas (|65ap and (|66p . 
we see that the numerator is a more sensitive function 
of the bias than the denominator, in particular diverging 
faster as 6 — >■ 0. We consequently divide our subsequent 
discussion into two cases: the strong bias regime in which 
(-'^tra) ^ (-'^fiuc), and the weak bias regime in which 
(-'i^tra) < (^fluc). We always assume (Xun^o) < (-'i^tra) 
in what follows, which can be insured by simply taking a 
sufficiently small choice of j]c in the definition of circular 
polarization at the beginning of Subsection IIV Al 
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a. Strong Bias Regime When (Xtra) ^ (^fluc), 
then the bias dominates the dynamics to the extent that 
the key distances characterizing the polarization dynam- 
ics are well-described by deterministic expressions. We 
will for the most part consider the typical case in this 
regime, in which (Xun^o) ^ (-'^fiuc) as well, and comment 
on what happens when this is not true later. Proceeding 
under the assumption that (Xun.o) 3> (^fluc), the stan- 
dard deviation of the distances Xnn.o and Xtra is com- 
parable to or smaller than their mean, so that these dis- 
tances are indeed comparable to the deterministic length 
scales Liin.o and Ltra with high probability. Moreover, 
the distance Xfluc is negligible relative to these other 
distances. Consequently, the switching region depth is 
approximately deterministic with Xdep ~ iiin.o + Lus.- 
This means that, with high probability, the soliton expe- 
riences exactly one polarization switch (as it moves into 
the favored polarization state), and does so in an approxi- 
mately deterministic manner after a distance Lun^Q+Ltm, 
as shown in the top panel of Fig. |51 When the ratio in 
Eq. (j67p becomes order unity, the polarization switch be- 
comes a bit more random (bottom panel of Fig. [5]), but 
usually multiple switches arc not seen. 

The case in which (Xnn.o) ^ (-^^fiuc) is only a minor 
modification of the above description, in that now Xh^^q 
behaves randomly but plays a negligible role in the dy- 
namics since necessarily {Xnnfi} ^ (-'^tra) given the def- 
inition of the strong bias regime. 

b. Weak Bias Regime For weak bias, when 
(-'^fiuc) ^ (-^tra), the distances Xna^ and Xtra become 
highly variable. This is to be expected since the limit of 
no bias 6 = involves a qualitatively different scenario 
described in Sec. IIVB II The weak bias regime involves 
features of both the no bias and strong bias regimes. On 
the one hand, the polarization will eventually collapse 
into a permanent circular polarization with the same 
sign as the bias b. However, for the case of weak bias 
we expect several or even many visits to the linear 
polarization state before the ultimate collapse into a 
circular polarization state. Indeed, in the limit of no 
bias, the linear polarization state is visited infinitely 
often, as discussed in Sec. lIVBip . 

The dynamics are in fact dominated by the extended 
period of random polarization switching since the initial 
approach toward the favored circular polarization state is 
short by comparison ((Xfluc) > (^iin,o>, (-'^tra)), though 
as noted above this initial approach has a highly random 
character because the standard deviations of Xiin,o and 
Xtra are large compared to their means. The polarization 
switching depth is therefore determined predominantly 
by the length scale Xjcp ^ Xa^c which has statistics 
described in Eq. (|SS|. Figure IHl illustrates how, as the 
bias is weakened (proceeding downward), the pulse po- 
larization switches randomly over an extended random 
distance Xfluc before finally collapsing into the favored 
polarization state. 
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FIG. 8. Realizations of the soliton polarization for a medium 
with strong bias (a = 0.1, b = —0.5, Lun.o = 0.64, I/tra = 4.26, 
ifluc ~ 0.04) on the top and a medium with not so strong bias 
(a = 0.3, b = -0.5, iiin,o = 0.64, Ltra = 4.26, Lruc = 0.36) 
on the bottom. 



V. CONCLUSIONS 

We have analyzed the resonant interaction of single- 
soliton light pulses with a A-configuration degenerate op- 
tical medium in the idealized integrable Maxwell-Bloch 
approximation. This is an example of a phenomenon for 
which integrability and structural disorder produce non- 
trivial stochastic nonlinear dynamics, yet whose statistics 
can be analyzed in closed form. We have found explicit 
dependence of the soliton polarization on the average dif- 
ference between the initial populations of the degenerate 
lower sub-levels along the medium sample, with infre- 
quent but persistent random switching between the two 
circular polarizations when this difference vanishes, and 
almost certain asymptotic approach to one of the two cir- 
cular polarizations determined by this average difference 
when it does not vanish. Moreover, we have provided a 
precise quantification of the statistical dynamics of po- 
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FIG. 9. Realizations of the soliton polarization for a medium 
with bias becoming progressively weaker from the top to bot- 
tom panel. Top: a = 0.5, b = —0.5, I/un.o = 0.64, Ltra ~ 4.26, 
Lfluc = 0.04; middle: a = 1.0, b = -0.5, Lun.o = 0.64, 
Ltra = 4.26, Lfluc = 4.00; bottom: a = 3.0, 6 = -0.5, 
Lun.o = 0.64, Ltra = 4.26, Lfl.c = 36.0. 



larization switching, including probability distributions 
for the key distances describing transitions. 

At least one question still remains about these re- 
sults, which is whether they are robust under the ran- 
dom fluctuations of the medium polarization induced by 
finitc-tcmpcrature effects. Our preliminary numerical re- 
sults confirm that they should be robust. These results. 



and also an analysis based on the full evolution equa- 
tions for the spectral data corresponding to random, non- 
vanishing initial medium polarization variables, as pre- 
sented in [48[, will be described in a subsequent publica- 
tion. 
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Appendix A: Correlation Length in the Optical 
Medium 



In this appendix, we give a mathematically precise de- 
scription of the correlation length, which is assumed to 
be effectively zero in the white noise approximation (|24p 
in Section fill Al In general, without the white noise as- 
sumption (j24bp . the correlation function of the popula- 
tion density difference a{x) in the medium is defined as 

R^{x) = {[aix') - b][a{x' + x) ~ b]), (Al) 

where we use the statistical spatial homogeneity of a{x). 
We further assume 



/"OO 

< / Raix) dx < OO, 
Jo 



(A2) 



which means that the correlations are sufficiently short- 
ranged. 

Using Eqs. (|A1[) and (jA2[) . the correlation length L^ 
is defined as 



1 r°° 1 f°° 

^c = TTTTTT / Ra{x)dx^^r Ra{x)dx. 

Ra(S>) Jo ^i Jo 



(A3) 



Physically, Lc is the shortest length scale over which the 
intial population density difference a{x) exhibits signifi- 
cant variations. 

Under the condition (IA2I) and for x ^ Lr,, the func- 



tional central limit theorem for random fields 58[ implies 
that the function A(x) = /g a{x)dx, defined in (|T7l) . is 
statistically equivalent at large x to the random process 
aW{x) -I- 6x, where h = {a{x)) is defined in p4a|) . 



Ra {x) dx 



1/2 



(A4) 



and W{x) is the standard Wiener process. 
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Appendix B: Probability distributions for transition 

SAvitching and interswitch distance in the absence of 

medium bias 



From the definition of Xtra, the formula (j23bp for the 
evolution of the ellipticity angle, and the white noise ap- 
proximation for the case of no medium bias {A{x) ~ 
aW{x)), we deduce 



X 



tra = min {x - xo : |?7(a;)| > 7r/4 - rjc} 

x>xo 



mm -^ X — Xo 

X>Xf) 



min \ X > Xo 

X>Xq 



^In^ 
2aT \d-\ 



W{x) 

- — r tanh~ (cos 2rir) > 

i\t\ '^'j 



2a 



m-) + :T^ln^ 



1 
2aT 



> V-^tra f- 



where in the above, xq is a position where ri{xo) = 
(equivalently W{xo) = -(l/2ar) ln(|d+|/|ci_|)). We see 
then that Xtra is just the distance of the first position 
X after xq at which the Wiener process W{x) escapes a 
given interval, given an initial position within that inter- 
val at Xo, also known as a first exit time. The formula 
(|57|1 then follows directly from the first exit time formula 
(2.8.24) in [6l| and the translational invariance of the 
statistics of Wiener process increments. 

The random variable Xjnt has a two-step definition, 
which in mathematical terms can be translated as follows: 



''^int — ^tra + -^^ret 



(Bl) 



where the "return distance" is defined as the distance 
over which the soliton returns from a nearly circular po- 
larization state of either orientation (|?7| ~ 7r/4 — 7]c) to a 
linear state (ry = 0): 



Xj-ot ~ min {a; — Xi : 'q{x) — 0} 

a:>a:i 

= min Ix — Xi : W{x) ~ — - — In — — - 
x>xi I 2aT \d-\ 



with xi a position where rj^xi) = 
alently W{xi) = -(l/2aT) ln(|d+| 



±(7r/4 — 77c) (equiv- 
'\d_\)±y/L^^). Be- 



cause of the statistical reflection symmetry of the Wiener 
process, either sign of the ± for the conditions at xi will 
give the same result, which is the first position x after 



xi at which the Wiener process W{x) achieves a certain 
value, given that it was situated at a different value at cci. 
This is known as a first passage time, and we can apply 
formula (2.8.5) in |6l| . again using statistical translation 
invariance, to obtain the probability density function for 



PX„t {x) 



2'!TX^ 



^-L,,j2x fQj. ^ > Q_ 



Now, by the strong Markov property of the Wiener 
process, the summands Xtra and X^-^t in Eq. (|Bip are 
independent random variables, so the probability den- 
sity of Xjnt is the convolution of the probability densities 
of the summands. Equivalently, the moment generating 
function 

/•OO 

Mx.„.(s)^(e^^-)= / e-px.„,(x)dx 



is the product of the moment generating functions of the 
summands |60| . These moment generating functions can 
be computed through the Laplace transform identity: 



:e"'''/*dt = e-^^^ for 6 > 0, (B2) 



V2M- 



which can be derived either by calculus tricks [6J| or, 
more elegantly, by computing the probability density 
function for the first passage time of Brownian motion by 
the refiection principle (Section 2.6 in [6l[) and its mo- 
ment generating function by the optional stopping the- 
orem (Section 2.8 in [61|) and connecting these results. 
Applying this Laplace transform identity, we obtain: 



x.„.(.) 


= Mx..A 


•)Mx.„ 


(^) 








OO 
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e{4"+i)v 


'Ltr^V2s^~^/ 


Ltra\/2i 




OO 

= 2>;o- 


-(4n+2)V 


J-traV^ _ 


— OO 

-2>;o'. 


VitraV^ 



)i=0 



But then we can infer the probability density (|57bp for 
^int by again applying the Laplace transform identity 
(IB2I) in reverse to each summand. 
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